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ABSTRACT 

Approximately 1/4 — 1/2 of short duration Gamma- Ray Bursts (GRBs) are fol- 
lowed by variable X-ray emission lasting ^ 100 s with a fluence comparable or ex- 
ceeding that of the initial burst itself. The long duration and significant energy of 
this 'extended emission' (EE) poses a major challenge to the standard binary neutron 
star (NS) merger model. Metzger, Quataert & Thompson (2008) recently proposed 
that the EE is powered by the spin-down of a strongly magnetized neutron star (a 
millisecond proto-magnetar) , which either survives the NS-NS merger or is created 
by the accretion-induced collapse of a white dwarf. However, the effects of surround- 
ing material on the magnetar outflow have not yet been considered. Here we present 
time-dependent axisymmetric relativistic magnetohydrodynamic simulations of the in- 
teraction of the relativistic proto-magnetar wind with a surrounding 10^^ — 10"'^ Mq 
envelope, which represents material ejected during the merger; in the supernova fol- 
lowing AIC; or via outflows from the initial accretion disk. The collision between the 
relativistic magnetar wind and the expanding ejecta produces a termination shock and 
a magnetized nebula inside the ejecta. A strong toroidal magnetic field builds up in 
the nebula, which drives a bipolar jet out through the ejecta, similar to the magnetar 
model developed in the case of long duration GRBs. We quantify the 'break-out' time 
and opening angle of the jet 0j as a function of the wind energy flux E and ejecta mass 
Mej. We show that E and dj are inversely correlated, such that the beaming-corrected 
(isotropic) luminosity of the jet (and hence the observed EE) is primarily a function 
of Mej. Both variability arguments, and the lower limit on the power of magnetar 
outflows capable of producing bright emission, suggest that the true opening angle of 
the magnetar jet must be relatively large. The model thus predicts a class of events 
for which the EE is observable with no associated short GRB. These may appear as 
long-duration GRBs or X-Ray Flashes unaccompanied by a bright supernova and not 
solely associated with massive star formation, which may be detected by future all-sky 
X-ray survey missions. 

Key words: magnetic fields - MHD - stars: neutron - star: wind, outflows - gamma- 
rays: bursts - method: numerical 



1 INTRODUCTION 

Gamma-Ray Bursts (GRBs) are canonically divided into 
'long soft' (LGRB) and 'short hard' (SGRB) classes, based 
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on their bimodal distribution in duration and spectral hard- 
ness (Kouveliotou et al. 1993). This division is supported 
by studies of the host galaxies and environments of GRBs. 
LGRBs occur in late type galaxies with high specific star for- 
mation rates (e.g. Fruchter et al. 2006; Levesque et al. 2010) 
and are accompanied by core collapse supernovae (SNe) 
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(Galama et al. 1998; Delia Valle et al. 2006; Chornock et al. 
2010; Starling et al. 2010). In contrast, SGRBs occur in both 
elliptical (Bloom et al. 2006; Berger et al. 2005) and late type 
(Fox et al. 2005; Barthelmy et al. 2005) galaxies; they show 
larger average offsets from their host centres (e.g. Procliaska 
et al. 2006; Berger 2009; Fong, Berger & Fox 2010); and no 
evidence is found for a bright associated supernova in a few 
well-studied cases (e.g. Hjorth et al. 2005; Kann et al. 2008). 

The differences between long and short GRBs have mo- 
tivated the development of distinct, albeit related, models 
for their central engines. LGRBs originate when the core 
of a rotating massive star collapses to form either a black 
hole (BH) (MacFadycn & Woosley 1999) or a rapidly spin- 
ning, strongly magnetized neutron star (a millisecond proto- 
magnetar, e.g. Usov 1992; Wheeler et al. 2000; Thompson, 
Chang & QuatEiert 2004). SGRBs are instead commonly at- 
tributed to the inspiral and merger of neutron star (NS)-NS 
or NS-BH binaries (Lee & Ramirez- Ruiz 2007; Nakar 2007), 
although the accretion induced collapse (AIC) of a White 
Dwarf (WD) (Metzger, Quataert & Thompson 2008) or a 
NS (Dermer & Atoyan 2006) represent viable alternatives. 
In all SGRB models, the short ~ 0.1 — 1 s duration of the 
burst is related to the accretion timescalo of the compact 
^ 10"'^ — O.IMq torus that forms around the central NS 
or BH following the merger (e.g. Janka et al. 1999) or AIC 
event (Dessart et al. 2006). 

The standard LGRB/SGRB dichotomy has recently 
been challenged by several 'hybrid' events that conform 
to neither class (e.g. Zhang 2007; Bloom, Butler & Perley 
2008). GRB 060505 and 060614 are both long bursts based 
on their duration, yet neither shows evidence for a bright 
associated SN (Fynbo et al. 2006; Gehrels et al. 2006; Gal- 
Yam et al. 2006; Ofek et al. 2007). Although the short GRB 
050724 occurred in an elliptical host with no associated SN, 
it was followed by variable X-ray emission lasting ~ 140 s 
with a total fluencc ~ 3 times greater than the short GRB 
itself (Barthelmy et al. 2005). Similar to 060614, the light 
curve of GRB 080503 was characterized by a hard initial 
spike, followed (after a brief lull) by a 'hump' of X-ray emis- 
sion lasting ~ 100 s and carrying ~ 30 times the fluence of 
the initial spike (Perley et al. 2009). The rapid variability 
of this 'extended emission' (EE) strongly suggests that it 
results from ongoing central engine activity. 

All together approximately ~ 1/4 of Swift SGRBs^ 
are accompanied by extended X-ray emission lasting for 
~ 10 - 100 s with a fluence ^ that of the GRB itself (sec 
Norris & Gehrels 2008 and Perley et al. 2009 for a compila- 
tion of events). The hybrid nature and common properties 
of these events ('short GRB' -|- ~ 100 s X-ray tail) have 
motivated the introduction of a new subclass: Short GRBs 
with Extended Emission (SGRBEEs). It was moreover re- 
cently discovered that some SGRBs are followed by an X-ray 
'plateau' ending in a very sharp break (GRB 980515; Rowl- 
inson et al. 2010). Troja et al. 2008; Lyons et al. 2010) and 
is difficult to explain by circumstellar interaction alone. Al- 
though the connection of this event to SGRBEEs is unclear, 



^ When observational bias due to the effects of e.g. detection 
threshold are taken into account, the true fraction of SGRBEEs 
could be as high as ~ 50% (Norris, Gehrels & Scaxgle 2010). 



it nevertheless provides additional evidence that the central 
engine is active at late times. 

An important observational question is whether SGR- 
BEEs differ from 'normal' short bursts in other properties. 
Troja et al. (2008) found that SGRBEEs occur on average 
closer to the centres of their host gala:xies than other SGRBs; 
however, the current paucity of well-localized events make 
statistical claims uncertain and recent studies have not ver- 
ified this result (Berger 2009; Fong, Berger & Fox 2010). 
Norris, Gehrels & Scarglc (2011) showed that SGRBEEs dif- 
fer also in the properties of the the initial short GRB itself. 
They find that the average duration of the burst, and of 
individual pulse structures, arc longer for short GRBs with 
EEs, possibly suggesting the existence of different progeni- 
tor channels (Leibler & Berger 2010; Virgili et al. 2011) or 
a different circumburst environment. 

The long duration and high fluence of the extended 
emission of SGRBEEs poses a serious challenge to the NS 
merger scenario, because in this model both the prompt 
and extended emission are necessarily powered by black 
hole accretion. It is in particular difficult to understand 
how such a high accretion rate is maintained at very late 
times. Although accretion of the torus formed from the 
merger may power the short GRB itself (e.g. Janka et al. 
1999), the timescale for the disk to accrete is short (typically 
< 1 s), and the disk is disrupted by outflows soon there- 
after (Metzger, Piro, & Quataert 2009a; Metzger, Piro, & 
Quataert 2009b; Lee, Ramirez-Ruiz & Lopez-Camara 2009). 
'Fall back' accretion could in principle power late emission 
(Rosswog 2007; Faber et al. 2006; Chawla et al. 2010), but 
whether sufficient mass is placed onto the highly eccentric 
(yet bound) orbits required by this scenario is unclear. Full 
simulations of the fall-back process, including the important 
effect of heating due to r-process nucleosynthesis (Metzger, 
Arcones, Quataert, & Martfnez-Pinedo 2010a), have yet to 
be performed. Even if sufficient mass returns at late times, 
the accretion will occur under radiatively inefficient condi- 
tions (e.g. Narayan, Piran & Kumar 2001); since the disk is 
only marginally bound when it cannot cool, most of the ac- 
creting mass may again be lost to [non-relativistic] outflows 
(Rossi & Begelman 2009). 



1.1 The Proto-Magnetar Model for SGRBEEs 

The hybrid properties of SGRBEEs hint that the engine 
could itself be of hybrid nature.^ Metzger, Quataert & 
Thompson (2008) recently proposed that SGRBEEs result 
from the birth of a rapidly spinning proto-magnetar, cre- 
ated by a NS-NS merger or the AIC of a WD (see Figure 1 
for an schematic illustration). In this model the short GRB 
is powered by the accretion of the initial torus (similar to 
standard NS merger models), but the EE is powered by a 



^ Lazzati, Morsoiiy & Bcgclman (2010) recently proposed that 
SGRBEEs may result from coUapsars viewed off-axis. However, 
the differences between the host galaxies and environments be- 
tween SGBEEs and long GRBs (Berger 2009; Fong, Berger & 
Fox 2010; Leibler & Berger 2010), and the lack in some events 
of an associated supernovae, indicates that this is not the main 
channel. 
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Figure 1. Schematic diagram of the stages of the proto-magnetar 
model for Short GRBs with Extended Emission. (A) The merger 
of two binary neutron stars, or the accretion-induced collapse 
of a rotating white dwarf, results in the formation of a compact 
^ 10^'^ — O.IMq torus around the central proto- neutron star. (B) 
Accretion of the torus powers a relativistic bipolar jet, resulting 
in a short GRB lasting ~ 0.1 — 1 s, similar to the standard NS-NS 
merger model. Following a<;cretion, however, a rapidly spinning 
(millisecond) proto-magnetar remains. (C) Material ejected dur- 
ing the merger, by the supernova following AIC, or via outflows 



from the ax;cretion disk, results in a 
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lope around the proto-magnetar moving outwards with a velocity 
Vf.j ~ 0.1 — 0.2 c. The relativistic wind from the proto-magnetar 
collides with the ejecta, producing a magnetar wind nebula. (D) 
Magnetic stresses in the nebula redirect the magnetar wind into a 
bipolar jet. After the jet breaks through the ejecta on a timescale 
~ 1 — 10 s (Fig. 3), the magnetar wind escapes and accelerates 
to ultraxelativistic speeds (Fig. 2). Emission from the jet at much 
larger radii powers the extended emission lasting ~ 10 — 100 s, 
similar to the proto-magnetar model for long GRBS (see Fig. 5). 



relativistic wind from the proto-magnetar at later times, af- 
ter the disk is disrupted. Although a NS remnant is guar- 
anteed in the CEise of AIC, the merger of a double NS bi- 
nary could also leave a stable NS remnant, provided that 
either (1) the total mass of the binary is low and/or the 
NS equation of state is stiff (Shibata & Taniguchi 2006); 
(2) the proto-NS forms in a meta-stable state supported by 
differential rotation (Baumgarte, Shapiro & Shibata 2000; 
Baiotti, Giacomazzo & Rezzolla 2008), but it then loses suf- 
ficient mass via magneto-centrifugal outflows (Thompson, 
Chang & Quataert 2004; Metzger, Thompson, & Quataert 
2007) to reach stability. The likelihood of this possibility has 
increased recently due to the discovery of a « 2Mq NS (De- 
morest et al. 2010), which suggests that the nuclear EOS is 
indeed stiff (sec also Ozol ct al. 2010). Given that rapid ro- 
tation is expected in both NS-NS merger and AIC scenarios, 
it is plausible that the proto-NS will generate a magnetar- 
strength field by, for instance, an a — dynamo (Duncan & 
Thompson 1992), shear instabilities at the merger interface 
(Price & Rosswog 2006) , or the magneto-rotational instabil- 
ity (MRI; e.g. Akiyama et al. 2003; Thompson, Quataert & 
Burrows 2005). 

Though not surrounded by the envelope of a massive 



star, magnctars formed from NS-NS mergers or AIC do 
not form in vacuum. In the AIC case ~ 10~^ — 10~^Mq 
is ejected during the SN explosion on a timescale < 1 s 
(e.g. Wooslcy & Baron 1992; Dessart et al. 2006), while in 
NS-NS mergers a similar mass may bo ejected dynamically 
due to tidal forces during the merger process (e.g. Rosswog 
2007). Mass loss also occurs in outflows from the accretion 
disk on timescales <^ seconds, due to heating from neutri- 
nos (Metzger, Thompson, & Quataert 2008; Dessart et al. 
2009), turbulent viscosity (Metzger, Piro, & Quataert 2008; 
Metzger, Piro, & Quataert 2009a), and nuclear energy re- 
leased by the recombination of free nuclei into ''He (Lee & 
Ramircz-Ruiz 2007; Metzger, Piro, & Quataert 2008; Lee, 
Ramirez- Ruiz & Lopez-Camara 2009). During the first few 
seconds after forming, outfiows from the magnetar itself axe 
heavily mass-loaded and non-relativistic, resulting in a sig- 
nificant quantity of ejecta ^ 1O~^M0 (Thompson, Chang & 
Quataert 2004; Bucciantini et al. 2006; Metzger, Thompson 
& Quataert 2007). All together, ~ 10~^ - Q.IMq is ejected 
with a characteristic velocity Uoj ~ 0.1 — 0.2 c and kinetic 
energy ~ 2 x 10^°(?;ej/0.1c)^(Mej/0.01MQ) ergs. 

A few seconds after the merger or AIC, one is left with 
a proto-magnetar embedded in a confining envelope.'* This 
configuration is qualitatively similar to that developed in 
the proto-magnetar model for LGRBs by Bucciantini et al. 
(2007, 2008, 2009), except that the enshrouding envelope is 
much less massive. In these previous works it was shown 
that, although the power in the magnetar wind is relatively 
isotropic (e.g. Bucciantini et al. 2006), its collision with 
the slowly-expanding ejecta produces a hot 'proto-magnetar 
nebula' (Bucciantini et al. 2007). As toroidal flux accumu- 
lates in the nebula, magnetic forces - and the anisotropic 
thermal pressure they induce - redirect the equatorial out- 
flow towards the poles (Begelman & Li 1992; Konigl & Gra- 
not 2002; Uzdensky & MacFadyen 2007; Bucciantini et al. 
2007, 2008, 2009; Komissarov & Barkov 2007). Stellar con- 
finement thus produces a mildly- relativistic jet, which drills 
a bipolar cavity through the ejecta. Once the jet 'breaks 
out', an ultra-relativistic jet (fed by the magnetar wind at 
small radii) freely escapes. The EE is then powered as the jet 
dissipates its energy at much larger radii. One virtue of ap- 
plying this picture to SGRBEEs is that it naturally explains 
why the EE resembles long GRBs in several properties, such 
as its duration and the existence of a late-time 'steep decay' 
phase (cf. TagHaferri et al. 2005; Perley et al. 2009). 

Although SGRBEEs resemble long GRBs in many prop- 
erties, important differences also exist. The EE is gener- 
ally softer (X-rays rather than gamma-rays) , somewhat dim- 
mer, and its variability is generally smoother (appearing to 
display e.g. a higher 'duty cycle'), than long GRBs. As- 
sessing the viability of the proto-magetar model for SGR- 
BEEs therefore requires determining whether these differ- 
ences may in part result from differences in the geometry 
of the relativistic outflow. These in turn may result because 



^ In cases when the ejecta originates from the earlier [non- 
relativistic] stage of the magnetar wind, the distinction between 
'wind' and 'ejecta' is blurred. In general, however, the magnetar 
outflow becomes ultra-relativistic relatively abruptly, such that 
this distinction is well-defined (Metzger et al. 2011). 



4 A^. Bucciantini, et al. 



the confining ejecta is significantly less massive and dense 
than in the core collapse case. 

In this paper we investigate the interaction of the rela^ 
tivistic proto-magnetar wind with the expanding ejecta us- 
ing axisymmetric (2D) rclativistic MHD simulations. Wc fo- 
cus in particular on the confining role of the ejecta and its 
dependence on the wind power, and on the ejecta mass and 
density profile. We show that collimation (jet formation) is 
achieved only within a bounded range of parameters. If the 
wind is too energetic, or the mass of the shell is too low, the 
ejecta is disrupted and little collimation occurs. In contrast, 
if the ejecta is sufficient massive and/or the wind is suffi- 
cient weak, the result is instead a 'choked jet' that may not 
emerge at all. We describe the numerical set-up in Section 
2 and present our results in Section 3. We apply our results 
to SGRBEEs in Section 4 and conclude in Section 5. 



2 NUMERICAL SETUP 

All calculations were performed using the shock-capturing 
central-scheme for rclativistic ideal MHD ECHO (Del 
Zanna, Bucciantini & Londrillo 2003; Del Zanna et al. 2007), 
using an ideal gas equation of state with an adiabatic coeffi- 
cient r = 4/3, as appropriate for relativistically hot gas. We 
refer the reader to these papers for a detailed description of 
the equations and numerical algorithms. 

We investigate the interaction of the magnetar wind 
with the surrounding ejecta envelope using 2D axisymmet- 
ric simulations on a spherical grid. The angular domain is 
9 = [0, it] with reflecting boundary at the polar axis to en- 
force axisymmetry, while the radial domain extends over the 
range r = [10^10^^1 cm. The grid in the radial direction 
is spaced logarithmically with 100 cells per decade, while 
spacing is uniform in the angular direction with 200 cells 
[we repeated selected simulations with twice the resolution 
to verify convergence; see also Camus et al. (2009) for es- 
timates of convergence with grid resolution, in similar sim- 
ulations as applied to pulsar wind nebulae (PWNe)]. We 
assume zeroth-order extrapolation at the outer boundary. 
The code is second-order in both space and time, with a 
monotonized central limiter, chosen in order to resolve the 
large density jump between the lighter relativistic plasma 
inside the magnetar wind nebula (MWN) and the heavier 
envelope. 

Bucciantini et al. (2007) showed that the interaction 
between a magnetar wind and a confining envelope depends 
on the strength of the toroidal magnetic field B in the 
MWN; this in turn depends on the magnetization in the 
wind (cr = r^B^c/E, where E is the wind energy flux) at 
the distance of the termination shock. While a can be cal- 
culated at the light cylinder radius with some confidence 
(e.g. Metzger et al. 2010), its value at larger radii is dif- 
ficult to determine due to uncertainties in the conversion 
of magnetic energy into kinetic energy in relativistic winds 
(see Bucciantini et al. 2007 for a detailed discussion of this 
problem in the context of LGRBs). Moreover, instabilities 
may occur inside the MWN (Begelman 1998) which further 
reduce the toroidal magnetic field strength. Nevertheless, 
previous studies (Bucciantini et al. 2008, 2009) show that 
reliable results for the dynamics of the MWN and the prop- 



erties of the jet are obtained using even the simplified regime 
of transverse relativistic MHD. 

For proto-magnetars with millisecond rotation periods, 
the light cylinder is located at ~ 10^ cm and the fast magne- 
tosonic surface is at ~ 10^ — 10** cm (Bucciantini ct al. 2006). 
As in the calculations of Bucciantini et al. (2008), at the in- 
ner boundary we inject a super-magnetosonic wind with a 
fixed Lorentz factor 7 = 25 and a magnetization of cr = 0.1; 
these values are appropriate for distances of the order of the 
termination shock radius if the conversion of magnetic to 
kinetic energy is efficient, and for typical properties of the 
magnetar wind at a few seconds after the proto-neutron star 
forms (Mctzgcr ct al. 2011; sec also the discussion in Sec. 3 
below). Under these assumptions, a is conserved throughout 
the upstream region. We assume that the wind contains a 
purely toroidal field and is cold with pc^/p = 100, where p 
and p are the density and pressure, respectively. For simplic- 
ity we assume that E and 7 arc constant, and that the wind 
is isotropic, throughout the simulation. Although in reality 
7 increases from ~ 1 to ~ 100 — 1000 over tens of seconds 
as the proto-magnetar cools (Metzger et al. 2011), studies 
of PWNe show that the dynamics becomes independent of 
7 in the limit 7^1. 

Following Darbha et al. (2010) wc assume that the shell 
of ejecta expands radially homologously with a velocity Vej = 
0.2c(r/re) inside a low density cavity of radius nn, where re 
is the radius of the outer edge of the shell. To this self-similar 
profile we further allow for a modulation with polar angle 6: 

_ j K{1 - acos2e)r~'^ if nn < r < re, 
\lO"^g cm-3 if r < nn. 

where re = 6 x 10^ cm, nn = 1.5 x lO^cm (corresponding 
to initial conditions one second after ejection), C = A, K 
is fixed by the total ejecta mass Moj, and we have cho- 
sen the density in the inner cavity to be sufficiently low 
as to have a negligible effect on the dynamics. The angu- 
lar profile we adopt is equatorially concentrated, motivated 
by the possibilities that either (1) the ejecta is the result 
of an cquatorially-focuscd disk wind (e.g. Metzger, Piro, & 
Quataert 2009a); or (2) a bipolar asymmetry remains from 
the jet produced during the early accretion-powered (short 
ORB) phase. Note that the self-similar profile of the ejecta 
implies that the ratio between the kinetic energy and mass 
of the ejecta is fixed at the value Eej/MejC^ — 0.005. Out- 
side re we Eissume a stationary, uniform medium with a den- 
sity p = 10~^g cm~^, which we have verified is sufficiently 
tenuous that it has negligible effects on the dynamics and 
geometry of jet formation. 



3 RESULTS 

Due to the self-similar nature of the ejecta, the dynamics 
of the wind-cjecta interaction is only a function of the ratio 
E/Mcj and the wind magnetization a, except as may result 
from latitudinal variations in the wind power. Even latitu- 
dinal effects are, however, suppressed by the formation of 
a MWN downstream of the termination shock. In this re- 
gion the wind is decelerated, plasma is heated, the magnetic 
field is compressed, and the pressure settles into a quasi- 
hydrostatic equilibrium, which depends only on cr. 

In Figure 2 we show the velocity and density profiles 
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Figure 2. Snapshots of the interaction of the proto-magnetar wind with the confining shell of ejecta at t = 35 seconds. The upper panel 
shows the density structure (cgs units), while the lower panel shows the magnitude of the velocity in units of the speed of light. The left, 
centre, and right columns show, respectively, cases corresponding to a low (E/M^j = 10^" erg ^p)^t case A), 'average' (E/M^j = 10^^ 



erg s 



M" 



case B) and high (E/M^j = 10^^ erg s ^ ^^Q^i case C) power wind. Axis are in units of lO'^^ cm. 



at t = 35 seconds from calculations performed for three dif- 
ferent values of the ratio E/Moj = 10^°, 10''\and 10^^ erg 
M!©^ (cases A, B, and C, respectively). Significant dif- 
ferences are immediately apparent in the properties of the 
outflow and the overall dynamics of the MWN-ejecta inter- 
action between the three cases. In both the low power (A) 
and 'average' (B) cases the MWN is confined within the 
ejecta. A well-coUimated bipolar relativistic outflow devel- 
ops, qualitatively similar to that found in the core collapse 
context as applied to LGRB (e.g. Bucciantini et al. 2008). In 
contrast, in the high power case (C) the MWN has almost 
completely blown the shell apart. 

Quantifying the geometry of the jet, in order to measure 
e.g. the jet opening angle, is nontrivial because its shape is 
not simply conical. In case A the jet shape is parabolic, while 
in case C the jet 'flares out' into a diverging flow. In the high 
power case C, it is unclear whether the ejecta will provide 
any confinement at all. Although Ray leigh- Taylor instability 
is fundamentally a three-dimensional process, and in princi- 
ple axisymmetric simulations might fail to reproduce prop- 
erly its detailed growth and geometrical properties, 2D sim- 
ulations in the context of Pulsar Wind Nebulae (Jun 1998; 
Bucciantini et al. 2004) agree with observations in term of 
the size and average properties of the unstable mixing layer. 



This might be due partly to the presence of a strong toroidal 
field which can suppress the growth in the azimuthal direc- 
tion. We note, moreover, that the jet propagates through 
the ejecta approximately an order of magnitude faster than 
the time required for the shell to fragment. Thus, even in the 
case of an energetic wind, a coUimated outfiow may form ini- 
tially in the polar region (albeit with a wide opening angle). 
In principle limited confinement could be hence maintained 
for a short time ~ 10 — 30 s, before shell fragmentation com- 
pletes and the outfiow becomes more isotropic. 

Figure 3 shows the 'break-out' time and characteris- 
tic opening angle of the jet as a function of E/Mcj, cal- 
culated from several simulations including those shown in 
Figure 2. The break-out time scales approximatively as 
oc {E/Mcj)~^^^ ■ Although, given the self-similar nature of 
the ejecta, we expect that quantities should depend primar- 
ily on the ratio E/Mej, the precise functional dependence 
is non-trivial to derive. Although we find that the basic jet 
properties are relatively robust to our assumed value for C, 
more substantial changes could in principle result for dif- 
ferent values of E/Mj, and rin. Nevertheless, we do not 
expect large variations in the latter quantities, with respect 
to the fiducial values adopted in this paper. 

Extrapolating Figure 3 to low values E/Mcj ^ 10*^ erg 
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Figure 3. Properties of tiie relativistic jet. Tiie upper panel 
shows the opening angle of the jet at t ~ 35 s (Fig. 2) as a 
function of the ratio E/M^j. Diamonds represent cases for which 
the ejecta shell is uniform with angle, while triangles show cases 
assuming a lower polar densities (a = 3 — 5; see cq. [1]). The lower 
panel shows the time that the jet 'breaks out' from the shell of 
ejecta (symbols are the same as in the upper panel). The dashed 
line is a power-law with exponent 0.5. 

Mq^, wc find that the jet requires > 20 s to break out. 
This timescale is comparable to both the delay observed 
before the onset of the EE in SGRBEEs, and to the time 
that the proto-magnetar wind spends at its highest spin- 
down luminosity (Metzger et al. 2011). E/M^j lO"'' erg 
s^^ thus represents a reasonable threshold power be- 

low which the jet is 'choked' inside the ejecta, in which case 
GRB-like high energy emission would fail. Naively it may 
appear that a similar criterion should apply also to the case 
of LGRBs, for which the typical mass of the surrounding 
envelope ^ lOMg, would imply a limit on the wind power 
> 10^*^ erg s~^. Note, however, several important differences 
between the dynamics of the MWN in the LGRB and SGR- 
BEE contexts. In the SGRBEE case, the ejecta expands at a 
significant fraction of the speed of light, such that at t ~ 20 s 
the shell radius (~ 10^^ cm) is significantly greater than the 
~ 10^^° cm radius of the progenitor envelope in the core col- 
lapse case. Since the polar expansion of the MWN is driven 
by internal pressure, the MWN suffers larger adiabatic losses 
in the SGRBEE scenario. The true threshold for a choked jet 
in the LGRB case thus occurs at lower luminosities > 10*^ 
erg s~^, more consistent with previous estimates in the lit- 
erature (Matzner 2003). 

Figures 2 and 3 show that the opening angle of the jet 
0j increases with E/Mej. One reason for this dependence 
is that 9j is proportional to the ratio of the pressure scale 
height of the MWN in the cylindrical direction H (resulting 
from magnetic stresses perpendicular to the vertical axis) to 
the radius of the ejecta Tc. Since H is itself proportional to 
the radius of the MWN (Bucciantini et al. 2007), which is 
larger for more energetic jets, it follows that Oj also increases 
with E/Mcj. 

In addition to the calculations performed for isotropic 
ejecta (as in Figure 2), in Figure 3 we also show results for 
cases with a lower polar density (a = 3 — 5 in eq. [1]). We 
find that the jet break-out times and opening angles show 
a similar dependence to the isotropic case, jjrovidcd that 
our results are parameterized in terms of an effective ejecta 
mass, defined as that of a spherically symmetric shell with 
a density equal to the polar value. 



We now briefly discuss the acceleration and variability 
of the jet. In all cases we find that the jet accelerates approx- 
imately linearly with radius (i.e. bulk Lorentz factor 7 oc r), 
suggesting efficient acceleration, as in traditional GRB fire- 
ball models (e.g. Goodman 1986; Paczynski 1986), although 
acceleration is slightly faster in the lower power case (A). 
This is consistent with previous studies of relativistic jets, 
for which it is found that acceleration is more effective for 
more narrowly collimated outflows (Komissarov 2010). We 
also find that the Lorentz factor of the jet has a 'top-hat' 
angular profile, with a constant-7 core and a sharp velocity 
gradient at the boundary to 7 > 1. 

Although the properties of the injected magnetar wind 
arc held constant in time, the properties of the jet as it 
leaves the ejecta are nevertheless variable. In particular, the 
Lorentz factor experiences order-unity fluctuations in the 
less energetic (more collimated) cases A and B on a length- 
scale of ~ 2 X 10^^ cm, corresponding to a typical timescale 
~ 5 — 6 sec. In the high power case C the jet properties show 
a smoother time evolution. One reason for this difference 
is that variability results from interaction in the magnetar 
nebula (in particular with the walls of the ejecta channel; 
Morsony, Lazzati & Bcgclman 2010): higher power jets are 
less variable because the ejecta channel is larger and pro- 
vides a less-coUimating environment. Indeed, we flnd that 
the observed timescales are similar to the sound/ Alfvenic 
time across the jet. Although we cannot rule out shorter 
timescale variability, as it would not be resolved, significant 
changes in the fluid geometry on small scales seem unlikely 
given the relatively large size of the MWN system. As we 
discuss in Section 4, because narrow jets axe predicted to be 
more variable than wider ones, the measured large amplitude 
variability of the EE provides an independent diagnostic of 
the jet opening angle. 

Because the proto-magnetar model for GRBs predicts 
a positive correlation between E and the wind magnetiza- 
tion a (Thompson, Chang & Quataert 2004; Metzger et al. 
2011), we also repeat the case B calculation for a higher 
value a = 0.2. We flnd that the jet break-out time is similar 
between the low and high magnetization cases, but the open- 
ing angle is a few degrees smaller in the high-a case. We want 
to stress here that we assume efflcint conversion of magnetic 
to kinetic energy either in the wind or at the termination 
shock such that the system behaves as if cr < 1 even for out- 
flows that are Poyting flux dominated at the Ligth Cylinder 
[see the discussion in Bucciantini et al. (2007) for the case 
of a proto-neutron star wind in the context of LGRBs]. The 
problem of magnetic dissipation in winds has been investi- 
gated in the case of pulsars (Kirk & Skjaeraasen 2003), while 
for for dissipation at the termination shock the reader is re- 
ferred to the recent results by Siroru & Spitkovsky (2011). 
Was also repeated our case B calculation for a shallower den- 
sity proflle (7 = 2 (eq. [1]). Again we found no signiflcant 
difference with respect to the reference case C = 4, which 
suggests that our results do not depend sensitively on the 
precise structure of the ejecta. 

4 IMPLICATIONS FOR SGRBEES 

After a clean opening is created through the ejecta, the 
power of the jet reflects, in a time- and angle-averaged sense. 
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Figure 4. Isotropic jet luminosity Ei^o = ^.f\^^> where = 
il/4n Si is tlie jet beaming fraction and Q is tlie opening solid 
angle of the jet, after it has relaxed on a tiinescale t ~ 35 s (Fig. 2). 
Diamonds show calculations performed assuming uniform ejecta, 
while triangles show cases with a lower polar density. 

the value of E{t) injected by the proto-magnetar wind at 
much smaller radii (e.g. Buceiantini ct al. 2009; Morsony, 
Lazzati & Begclman 2010). Most emission models for GRBs, 
such as internal shocks or magnetic reconnection, are the 
result of internal dissipation in the jet. In this case the ob- 
served (photon) luminosity is typically proportional to the 
beaming-corrected, or isotropic, jet luminosity £'iso = Ef^^, 
where /b ^ /2 is the beaming fraction. 

Figure 4 shows Eiao as a function of Mej / E for the same 
calculations shown in Figure 3. The inverse correlation be- 
tween the wind power and opening angle (Fig. 3) implies 
that, at fixed ejecta mass, Eiso is constant to within a fac- 
tor ~ 3 across two orders of magnitude in E. Thus, as long 
as the shell effectively coUimates the outflow, the EE lumi- 
nosity in the proto-magnetar model depends primarily on 
Afej. If, by contrast, E/Mcj > W^'^ erg s"^ Mq^ then the 
shell may be entirely disrupted (in which case the isotropic 
luminosity is instead directly oc E), while if E/Mgj < Iff*® 
erg s~^ Mq* the jet is probably choked and no emission is 
expected on timcscales of relevance. 

The geometry of the magnetar jet also has consequences 
for the ubiquity of EE associated with short GRBs. As noted 
by Metzger, Quataert & Thompson (2008), without con- 
finement the magnetar outflow (responsible for the EE) is 
mostly equatorial, while the accretion-powered jet (respon- 
sible for the initial short GRB) is probably polar. An im- 
portant question is thus whether a typical observer will see 
both components. Our results show that, except perhaps 
in the most energetic cases, the magnetar wind is diverted 
into a polar outflow. Unfortunately, the opening angles of 
SGRBs are poorly constrained* observationally, with mea- 
sured values ranging from a few to > 25 degrees (Burrows 
ct al. 2006; Grupc et al. 2006). Recent numerical results 
suggest values in the range of ~ 20 degree (RezzoUa et al. 
2011). It is thus possible that events could exist for which 
the extended emission is not observable because it is more 
coUimated than the initial SORB, in which case the event 
would be classified as a 'normal' short burst (see Barkov & 

* In such a 'two jet' scenario it is also unclear which jet to as- 
sociate a putative opening angle measurement with (e.g. Granot 
2005). 



Pozanenko 2011 for a similar idea). Such events carmot be 
too common because the fraction of short GRBs with ob- 
served EE is already rather large (Norris &: Bonnell 2006). 
This implies that the magnetar wind cannot be too coUi- 
mating, which suggests that the average shell mass is low 
(we provide additional evidence for a wide-angle magnetar 
jet below). 

Alternatively, events may exist for which only the EE 
is observable, because the initial short burst is more nar- 
rowly coUimatcd. These events would probably be classified 
as regular long duration GRBs or X-ray Flashes, but would 
not be accompanied by a bright associated supernova. It is 
difficult to place definitive constaints on the rate of such 
events, although we note that at least one X-ray Flash with 
an EE-like hght curve was not in fact accompanied by a 
bright supernova (XRF 040701; Soderberg et al. 2005). 

We now attempt to constrain the properties of the 
ejecta using the measured luminosity of the EE. The sam- 
ple of SGRBEEs with known redshifts and measured EE 
fluences is unfortunately small and incomplete. The sample 
may furthermore be biased against less luminous events, in 
which case the lower limits are not constraining. Neverthe- 
less, when measured, the isotropic luminosity of the EE is 
typically in the range Lee ~ 2 x 10"'** — 2 x 10"'® erg s~^ 
(Figure 5 shows some examples). Since we found that dur- 
ing the early jet-formation phase the isotropic luminosity is 
(Fig. 4) 

Sisoj ~ 1 - 3 X lO^'(Mej/O.lM0)crg s"\ (2) 

we can relate the observed EE luminosity Lee to the ejecta 
mass: 

Mej ~ 0.01- 0.03 ("—j^^^") (Il^y (ni^y Mr., 
' V10*®ergs-V ^0.1/ V 0.3 / ^ 

(3) 

where rjrad = L-be/EsoME is the radiative efficiency of the 
jet, and r/p = iJiso.EE/Sisoj ~ 0.1 — 0.3 is the ratio between 
the isotropic power of the magnetar wind during the EE 
phase i^iso.BE at late times (t ~ 10 — 100 s) and that at 
early times Eiso,] {t ^ 10 s), when the opening angle of the 
jet is determined. Detailed evolutionary models of proto- 
magnetar spin-down (Metzger et al. 2011) show that r?p is 
typically ~ 0.1 - 0.3. 

Equation (3) shows that for typical values of r?rad, ^ee, 
and the measured EE luminosity, the inferred ejecta masses 
are in the range Mcj ~ 10~^ — 10~^Mq, consistent with the 
range of predicted ejecta masses in both the NS-NS merger 
and AIC scenarios (Sec. 1). This represents an important 
consistency check on the proto-magnetar model. 

Adopting a typical value of the ejecta mass Mej = 
10~'^Mq, our results require that the early-phase power 
of the proto-magnetar wind must lie in the range E ~ 
10** — 10^° ergs s~^ so as to both maintain confinement and 
jet coUimation, yet not produce a choked jet. This spin-down 
luminosity corresponds to magnetars with either a relatively 
slow initial rotation rate (spin period Po ^ 3 ms) or a rela- 
tively weak dipole magnetic field strength B^ip < 2 x lO'^^ 
G (Metzger et al. 2010), as compared to those required to 
produce a bright Classical GRB. However, rotation periods 
much larger than the break-up rate (Po ~ 1 — 2 ms) seem 
unlikely given the substantial angular momentum that is re- 
quired in both NS-NS merger and AIC scenarios to form an 
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Figure 5. Average bolometric luminosity of GRB emission from 

the proto-magnctar jet as a function of time after formation, 
calculated using the models described in Metzger et al. (2011). 
Emission predicted by the internal shock and magnetic dissipa- 
tion models are shown with solid and dashed lines, respectively. 
The calculation assumes that the magnetar has an aligned dipole 
field of strength Bjip = 2 x lO^^ G and an initial spin period 
Pq = 1.5 ms. We adopt a value for the electron radiative efficiency 
Ee = 0.2 and a beaming fraction /;, = 0.3 (see text). For compar- 
ison we also plot the 15 — 350 keV Swift BAT extended emission 
light curves for GRBs 060614 (dotted), 080503 (dot-dashed), and 
061005 (triple-dot-dashed) (Butler & Kocevski 2007). 

initial disk. Thus, we are led to consider somewhat lower 
field magnotars. 

In Figure 5 we show models for the average bolomet- 
ric luminosity of GRB emission from the proto-magnetar 
jet, calculated using the models described in Metzger et al. 
(2011). We show the predicted emission in both internal 
shock and magnetic dissipation models, assuming the mag- 
netar has an aligned dipole field of strength Bdip = 2 x 10^^ 
G and an initial spin period Po = 1.5 ms. For comparison, 
we also plot the extended emission light curves of several 
Swift-detected SGRBEEs. In this model the jet power at 
f ~ 1 - 10 s is S ~ 10"'' - 10''" erg s"\ Applying these 
values to Figure 3, we adopt a beaming fraction /b ~ 0.3 
corresponding to a relatively large opening angle 6j ~ 45° . 

Figure 5 demonstrates that the emission predicted by 
the proto-magnetar model qualitatively agrees with the ob- 
served onset, duration, and luminosity of SGRBEEs. Also 
note the presence of a phase of steep decay in some of the 
SGRBEE light curves at late times. A relatively abrupt shut 
off in prompt emission is indeed expected to occur in the 
magnetar model when the (initially opaque) proto-neutron 
star becomes transparent to neutrino emission on a timescale 
~ 30 — 100 s (see Metzger et al. 2011 for a discussion). The 
observed spectral properties of SGRBEEs can also be com- 
pared to the predictions of our emission models; in GRB 
060614, for instance, the spectrum peaked at an energy 
i?poak ~ 50 kcV, which decreased by a factor ~ 2 throughout 
the EE (e.g. Zhang et al. 2007). For internal shock emission, 
the value of Speak and its time evolution are unfortunately 
very sensitive functions of the assumed microphysical pa- 
rameters, which makes a definitive comparison to our model 
challenging. By contrast, the magnetic dissipation model in- 
deed predicts a lower value Speak <, 100 keV for jets with 



smaller values of Eiso and a (see Metzger et al. 2011; their 
eq. [11]) as compared to those associated with standard long 
duration GRBs. 

Although Figure 5 shows a reasonably good fit to the 
EE light curves, similar agreement would not obtain for 
lower-i?(jip magnetars with less powerful winds. Less pow- 
erful outflows are 'dirtier' and require at least several tens 
of seconds before becoming optically-thin at the emission 
radius (see Fig. 8 of Metzger et al. 2011), inconsistent with 
the much earlier observed onset of the EE. Because the mag- 
netar wind cannot be too weak, self- consistency with the 
results of Figure 3 demands that the jet opening angle can- 
not be too small. We want to stress here that the problem 
of relating the outflow kinetic luminosity with the observed 
7-ray emission, depends on assumptions about acceleration 
efficiency and radiation mechanism that are still poorely un- 
destood, even in the canonical collpasar model of LGRBs. 
The present rcsuls should be considered indicative that the 
model can in principle reproduce the data at least in their 
broad temporal evolution. Indeed, the magnetar wind could 
also be much more powerful powerful (E 10*" ergs s~^), in 
which case the shell may be completely blown apart (case C 
in Fig. 2). As discussed above, in this case the rate of events 
for which the EE would be observable without an accompa- 
nying short GRB would be quite high. Future all sky X-ray 
missions (e.g. JANUS; Fox & JANUS Team 2010) could de- 
tect SGRBEEs by triggering on their EE alone, and hence 
may capable of addressing this question. 

An independent way to break the degeneracy between 
wide and narrow-angle jets is with variability. As discussed 
in Sect. 3, the amplitude of fluctuations in the jet proper- 
ties scale inversely with the jet opening angle. Since EE are 
generally less variable than standard long duration GRBs, 
and the latter are estimated to have opening angles in the 
range of a few-10 degrees (e.g. Frail et al. 2001), this also 
hints that larger opening angles are favoured for SGRBEEs. 



5 CONCLUSIONS 

The long durations and other similarities of the extended 
emission following some short GRBs is difficult to under- 
stand within the standard NS merger scenario. Building on 
the work of Metzger, Quataert & Thompson (2008), in this 
paper we present numerical simulations of the interaction 
between the energetic wind from a remnant proto-magnetar 
and the ejecta from the merger or AIC event. We analyze 
the confining properties of the ejecta shell and calculate the 
dependence of the jet properties on the properties of the 
proto-magnetar wind. 

Our results show that, for the self-similar ejecta struc- 
ture that we have assumed, the evolution of the system is 
primarily a function of the ratio between the power of the 
magnetar wind and the mass of the ejecta. Latitudinal den- 
sity variations can be reparametrized in term of an effective 
mass. We also demonstrate that well-collimated jets form 
only over a limited range in energy: a sufficiently powerful 
wind completely disrupts the confining shell, while a lower 
wind instead produces a choked jet. Interestingly, we find an 
anti-correlation between the wind power and opening angle 
(in the range of E/Mei over which a collimated jet forms), 
which implies that the observed (isotropic) luminosity is al- 



SGRBEEs from proto-magnetar spin-down 9 



most independent of E/Mcj (Fig. 4). Adopting a fiducial 
shell mass Mcj ~ 10^^ Mq, the corresponding luminosity 
~ 10"'^ ergs s~^ is broadly consistent with observed EE lu- 
minosities (Fig. 5). 

Considering the wind properties necessary to produce 
bright non-thermal emission also constrains the types of 
magnetars capable of powering the EE. Only magnetars with 
relatively strong surface dipole fields Saip ^ 10^^ G produce 
winds that are sufficiently 'clean', such that they become 
optically thin at the radii where the jet energy is dissipated 
sufficiently early to explain the ~ 10 s observed onset of the 
EE. Consistent with this notion, we find that a millisecond 
magnetar with > 10^^ G indeed produces emission qualita- 
tively consistent with observed EE light curves (Fig. 5). 

Our investigation has primarily focused on a fiducial 
profile for the ejecta. As discussed in the Introduction, al- 
though the ejecta mass may vary by a factor ^ 100 de- 
pending on e.g. the scenario responsible for producing the 
magnetar (Fig. 1), the characteristic size and velocity of the 
ejecta are not expected to vary by more than a factor of a 
few. For these reasons our results may be relatively robust. 
We caution, however, that a full parameter study is neces- 
sary to properly investigate questions such as to what degree 
the opening angle and break-out time depend on additional 
parameters, such as the relative size of the inner and outer 
radius of the ejecta. 

One limitation of our calculations is that the properties 
of the jet depend on the assumed magnetization of the wind 
a = 0.1, which was chosen because (1) previous results in 
the context of LGRBs (Bucciantini et al. 2008) agree with 
more detailed calculations where the magnetization in the 
wind evolves according to the value near the light cylin- 
der predicted from evolutionary models including PNS cool- 
ing (Bucciantini et al. 2009); and (2) because a variety of 
processes at work in the wind, the termination shock, and 
the nebula, may dissipate the toroidal flux, as is inferred 
in pulsar wind nebulae. Recently Mizuno et al. (2011) have 
shown that current driven instability can efficiently reduce 
the amount of toroidal magnetic field. They conclude that 
the magnetic field in the nebula can behave as if the eflfective 
a is much lower (less by a factor 10) than the true value. 
If this also holds in the more complex scenario of a MWN, 
then it is possible that current-driven instabilities may reg- 
ulate the magnetization, such that the lower effective value 
of cr as we have chosen is appropriate. 

Up to this point we have not discussed how to distin- 
guish between proto-magnetars formed from AIC and via 
NS-NS mergers, since the physical scenario is similar (see 
Fig. 1). Although in neither case do we expect the event 
to be associated with bright supernova, the ejecta is nev- 
ertheless composed of radioactive isotopes. As this material 
decays to stability, it reheats the ejecta, powering transient 
optical emission lasting ~ 1 day, with a peak brightness 
pLi, ^ ]^q41-42 ^j.gg approximately a thousand times 
brighter than a nova (Metzgcr, Piro & Quatacrt 2009b; Met- 
zger ct al. 2010b). Although such 'kilonova' have not yet 
been definitively detected (although see Perley et al. 2009), 
current searches are underway (Kasliwal et al. 2010). One 
way to distinguish between AIC and NS-NS mergers is that, 
in AIC the ejecta is predicted to be rich in Fe-group elements 
(e.g. ^®Ni; Metzger, Piro & Quataert 2009b; Darbha et al. 
2010), while in NS-NS mergers the ejecta is composed of 



more exotic, heavy r-process nuclei (Metzgcr et al. 2010b). 
We note that if the proto-magnetar injects additional en- 
ergy behind the expanding ejecta, and if this energy is ther- 
malized, this could enhance the luminosity of the predicted 
kilonova emission, as in analogous models for very luminous 
core collapse supernova (Kasen & Bildsten 2010; Woosley 
2010; Piro & Ott 2011). We plan to explore this possibility 
further in future work. Evidence of a large opening angle for 
magnetar outflows in SGRBEEs suggests the possibility of a 
strong interaction between the flow and (1) the companion 
in the AIC scenario, (2) a possible tertiary component in the 
NS-NS merger scenario (Thompson 2010), or (3) the clumps 
of matter ejected in unbound orbits in both scenarios. Such 
interactions could lead to observational signatures, which 
might help discriminate among the various scenarios (Mari- 
etta, Burrows & Fryxell 2000; MacFadyen, Ramirez-Ruiz & 
Zhang 2005). 

Another characteristic distinguishing NS-NS mergers 
and AIC is the strength and form of their gravitational wave 
signal (Abdikamalov et al. 2010). The coincident detection 
(or constraining upper limits) on gravitational wave emis- 
sion from SGRBEEs with future detectors such eis Advanced 
LIGO/Virgo should conclusively settle this issue. The detec- 
tion of a giant SGR-like flare from the location of a previ- 
ous SGRBEE (Giannios 2010) would be 'smoking gun' proof 
that magnetar birth is at the origin of these events. 
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